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1 Introduction: Bags and Resolvents 



A central concept in particle physics states that fundamental particles acquire their 
masses through interactions with vacuum condensates. Thus, a massive particle may 
carve out around itself a spherical region |]J or a shell |2| in which the condensate is 
suppressed, thus reducing the effective mass of the particle at the expense of volume 
and gradient energy associated with the condensate. This picture has interesting 
phenomenological consequences [[[], |3|. 

This phenomenon may be studied non-perturbatively in model field theories in 
1 + 1 space-time dimensions such as the Gross-Neveu (GN) model || and the multi- 
flavor Nambu-Jona-Lasinio (NJL) [f| model, in the large N limit. 

Explicit calculations of fermion bag profiles in the GN and NJL models were given 
originally in ||, [|7| and in ||. 

Following these works, fermion bags in the GN and NJL models were discussed in 



the literature several other times |§, flO ], ||L1 |, ]12] . For a recent review on these and 



related matters, see HT3 



Very recently, static chiral fermion bag solitons in a 1 + 1 dimensional model, 
as well as non-chiral (real scalar) fermion bag solitons fl5l , were discussed, in which 
the scalar field that couples to the fermions was dynamical already at the classical 
level. 

Mathematical considerations similar to those involved in studying fermion bags, 
appear also in other branches of theoretical physics, such as the theory of inhomoge- 
neous superconductors flEfl , and the results of this paper may be applicable there as 
well. 

Studying the physics of fermion bags necessarily involves knowledge of the re- 
solvent of the Dirac operator in the background of the bag. As an example, let us 
consider the 1 + 1 dimensional NJL model (which contains the GN model as a special 
case). 
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One version of writing the action of the 1 + 1 dimensional NJL model is 

S = Jd 2 x ij2 M ~ ( a + m ^)] % - 2^ + vr 2 )| , 

where the ip a (a = 1, . . . , N) are iV flavors of massless Dirac fermions, with Yukawa 
couplings to the scalar and pseudoscalar auxiliary fields a(x),it(x). 
The partition function associated with ( |1.1|) is 



Z = / Da Vtt Dip Vtp exp i / d 2 x I i/j 



{a + iir^) 



-^(^ 2 )} (L2) 



Integrating over the grassmannian variables leads to Z — J T>aT>n exp{iS e ff[a, 7r]} 
where the bare effective action is 

S eff [a,n] = -—J d 2 x (V + vr 2 ) - iN Tr log [i0 - (a + mj 5 )] (1.3) 

and the trace is taken over both functional and Dirac indices. 

This theory has been studied in the limit N — > oo with Ng 2 held fixed 0]. In this 
limit, ( |1.2j ) is governed by saddle points of (|1.3| ) and the small fluctuations around 
them. The most general saddle point condition reads 



SS eff 

5a (x, t) 



a (x, t) 



+ iN tr 



(x, t\ 



10- [a 



Z7T7 5 ) 



\x,t) 



5S 



eff 



71 (x, t) 



Stt (x, t) 



- N tr 



75 (x,t\ 



\x,t) 



0. 



- (a + Z7T7 5 ) 

Fermion bags are the space-time dependent solutions (a (x, t) , tt (x, t)) of ( |1.4j ), sub- 
jected to appropriate boundary conditions at spatial infinity, and on which S e ff/N is 
finite. 

Thus, studying fermion bags necessarily involves the resolvent of the Dirac oper- 
ator in the background of the bag. 

In this paper we discuss some mathematical aspects of the much simpler problem 
of static fermion bags, namely, the static solutions (<r(x), 7r(x)) of ( |OD . 
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For the usual physical reasons, we set boundary conditions on our static fields 
such that a(x) and ir(x) start from a point on the vacuum manifold a 2 + it 2 = m? 
(with constant a, tt of course, and where m is the dynamical mass M ) at x — — oo, 
wander around in the a — it plane, and then relax back to another point on the vacuum 
manifold at x = +00. Thus, we must have the asymptotic behavior 

a — > mcos^i , a 1 — > 

x^±oo x— >±oo 



71 — > msin#± , 7r' — > (1-5) 

x— >±oo x— >±oo 

where 9± are the asymptotic chiral alignment angles. Only the difference 9 + — 9- is 
meaningful, of course, and henceforth we use the axial U(l) symmetry of ( |1 . 1| ) to set 
9- = 0, such that a(— 00) = m and 7r(— 00) = 0. We also omit the subscript from 
9 + and denote it simply by 9 from now on. As typical of solitonic configurations, we 
expect, that cr(x) and it(x) tend to their asymptotic boundary values ( |1.5| ) on the 
vacuum manifold at an exponential rate which is determined, essentially, by the mass 
gap m of the model. 

Thus, in order to study static fermion bags, we need to invert the Dirac operator 

D = i$ - 00) + i7r(x)7 5 ) (1.6) 
in a given background of static field configurations <y{x) and tc(x), subjected to the 



boundary conditions (|1.5| ). In particular, we have to find the diagonal resolvent of 
( |1.6| ) in that background. We stress that inverting ( |1.6| ) has nothing to do with the 
large iV approximation, and consequently our results are valid for any value of N. 

The rest of the paper is organized as follows: In Section 2 we show that the 
Dirac equation (i0 — 00) + i7r(s)7 5 )^ ip = in a given static <j(x) + ry 5 7rO) back- 
ground, is equivalent to a pair of two isospectral Sturm-Liouville equations in one 
dimension, which generalize the well known one-dimensional super symmetric quan- 
tum mechanics. We use this generalized supersymmetry to express all four entries of 
the space-diagonal Dirac resolvent (i.e., the resolvent evaluated at coincident spatial 
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coordinates) in terms of a single function. In Section 3, we use the results of Section 
2 to derive simple expressions for various bilinear fermion condensates in the given 
static a(x) + ij 5 n(x) background. In particular, we prove that each frequency mode 
of the spatial current (ip(x)^ l ip(x)) vanishes identically in the static background. 
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2 Resolvent of the Dirac Operator With Static 
Background Fields 



As was explained in the introduction, we need to invert the Dirac operator ( |1.6| ). 
D = i$ — (cr(x) + i7r(x)7 5 ), in a given background of static field configurations a(x) 
and n(x), subjected to the boundary conditions (|I~5| ). 
In this paper we use the Majorana representation 

7° = 0"2 , 7 1 = ios and 7 5 = — 7°7 1 = <J\ (2.1) 

for 7 matrices. In this representation (|1.6|) becomes 



-d x — a —iuj — m \ I —Q —iu — in \ 
D= I = " , (2.2) 

iuj — Z7T d x — a J \ iuj — iir —Q^ J 

where we introduced the pair of adjoint operators 

Q = a(x) + d x , Q^ = a(x)-d x . (2.3) 



(To obtain ( p.2|) , we have naturally transformed %$— [p{x) + m{x)^ to the tu plane, 
since the background fields <j(x),tt(x) are static.) 
Inverting ( [2.2|) is achieved by solving 



/ —Q —iu — m{x) 
\ iuj — in(x) —Q^ 



I a(x,y) b(x,y) 
V c(x,y) d(x,y) 



-ilS(x-y) (2.4) 



for the Green's function of (|2.2|) in a given background a(x) } ir(x). By dimensional 



analysis, we see that the quantities a, b, c and d are dimensionless. 

2.1 Generalized "Supersymmetry" in a Chiral Bag Back- 
ground 



We now show that the spectral theory of the Dirac operator (|2.2j ) is underlined by a 
certain generalized one dimensional supersymmetric quantum mechanics. This gen- 
eralized supersymmetry is very helpful in simplifying various calculations involving 
the Dirac operator and its resolvent. 



The diagonal elements a(x,y), d(x,y) in ( |2.4j ) may be expressed in term of the 
off-diagonal elements as 



a(x,y) 



—i 



Q^c(x,y), d(x,y) 



uj — tt(x) lu + it(x) 

which in turn satisfy the second order partial differential equations 
f 



Qb(x,y) (2.5) 



Q 



Q 



UJ + 7r(x) 

d x b(x, y) 
uj + ir(x) 

1 



Q - {uj - 7r(x)) 



b(x,y) 



+ 



a 



UJ + 7r(x) 



uj + ir(x) 



5{x - y) 



-d x 



UJ — 7l(x) 

d x c(x, y) 



g f - (UJ + 7T{X)) 



c(x,y) 



UJ — 7r(x) 



+ 



a(x) 2 + vr(x) 2 + a'(x) - uj 2 + 



OKX 7T \X\ 



UJ — 7l(x) 



c(x,y) 

UJ — 7r(x) 



-5(x - y) . 
(2.6) 



Thus, b(x, y) and —c(x, y) are simply the Green's functions of the corresponding 
second order Sturm-Liouville operators 1 

b(x) 



L b (u)b(x) = -d x 



LJu)c(x) = -d x 



in (p.6|), namely, 



d x b(x) 


+ 


uj + ir(x) 


d x c(x) 


+ 


uj — ir(x) 



a(x) 2 + n(x) 2 - a'(x) -u 2 + 

UJ + Tl(x) 



a(x) 2 + tt(x) 2 + a'(x) - u 2 + a ( x W( x ) 



UJ — 7r(x) 



uj + ir(x) 
c(x) 

UJ — 7t(x) 

(2.7) 



b(x,y) 



9(x-y) b 2 {x)b 1 (y) + 9 (y - x) b 2 (y)bi(x) 
W h 



1 Note that u> plays here a dual role: in addition to its role as the spectral parameter (the w 2 
terms in (2.7)), it also appears as a parameter in the definition of these operators- hence the explicit uj 
dependence in our notations for these operators in (2.7). For this reason, it is possible to completely 



factorize the operators Lt and L c by additional obvious w-dependent similarity transformations on 
Q and QK However, these similarity transformations are singular at points where ir(x) — ±uj and 
are thus ill defined, and we will avoid them. 
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c ( x y j = 9 (g ~ V) c 2(x)ci(y) + 8(y-x) c 2 (y)c 1 (x) ^ g ^ 

Here {bi(x),b 2 (x)} and {ci(x), c 2 (x)} are pairs of independent fundamental solutions 
of the two equations L b b(x) = and L c c(x) = 0, subjected to the boundary conditions 

h(x), Cl {x) — 4 1 c(^)e _te > & 2 (x),c 2 (x) — A, c (A;)( 2 )e^ (2.9) 

X— > — oo ' rr^+oo 

with some possibly k dependent coefficients A^l(k), A^l(k) and withQ 



k = Vuj 2 - m 2 , lm£;>0. (2.10) 

The purpose of introducing the (yet unspecified) coefficients A^l(k), Af\(k) will be- 
come clear following Eqs. (|2.13|) and (|2.14|) . The boundary conditions ( |2.9|) are 
consistent, of course, with the asymptotic behavior ( |1.5|) of o and tt due to which 
both Lb and L c tend to a free particle hamiltonian [—d 2 + m 2 — u 2 } as x — > ±oo. 
The wronskians of these pairs of solutions are 

b 2 (x)b' 1 (x) — bi(x)b' 2 (x) 



W b (k) 
W c {k) 



UJ + 7r(x) 

c 2 (x)c[(x) - c 1 (x)c' 2 (x) 

UJ — 7r(x) 

(2.11) 



As is well known, W b {k) and W c (k) are independent of x. 

Note in passing that the canonical asymptotic behavior assumed in the scattering 
theory of the operators L b and L c corresponds to setting A bc = A bc = 1 in ( |2.9|) . 



Thus, the wronskians in (|2.11|) are not the canonical wronskians used in scattering 



theory. As is well known in the literature , the canonical wronskians are propor- 
tional (with a k independent coefficient) to k/t(k), where t(k) is the transmission 
amplitude of the corresponding operator L b or L c . Thus, on top of the well-known 



2 We see that if Im/c > 0, b\ and c\ decay exponentially to the left, and 62 and C2 decay to the 
right. Thus, if Xmk > 0, both b(x,y) and c{x,y) decay as \x — y\ tends to infinity. 
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features of t(k), such as the fact that t(k) has simple poles on the positive imaginary 
fc-axis (corresponding to bound states), the wronskians in (|2.11|) will have additional 
spurious fc-dependence coming from the amplitudes A^l(k), A^l(k) in (|2.9|). 

Substituting the expressions (|2.8|) for the off-diagonal entries b(x, y) and c(x, y) 
into ( |2.5| ), we obtain the appropriate expressions for the diagonal entries a(x,y) and 
d(x,y). We do not bother to write these expressions here. It is useful however to 
note, that despite the <9 x .'s in the Q operators in (|2.5|) , that act on the step functions 
in (|2~8|), neither a(x,y) nor d(x,y) contain pieces proportional to 5{x — y) . Such 
pieces cancel one another due to the symmetry of (|2.8|) under x <-> y. 

We will now prove that the spectra of the operators L b and L c are essentially the 
same. Our proof is based on the fact that we can factorize the eigenvalue equations 
L b b(x) = and L c c[x) = as 

1 Q f 1 -—Qb = b 



CO — 7r(x) CO + 7t(x) 



u) + ir(x) uj — 7r(x) 

(2.12) 

as should be clear from ( |2.6| ) and ( |2.7| ). 

The factorized equations ( |2. 12| ) suggest the following map between their solutions. 
Indeed, given that Lf,b(x) = 0, then clearly 

c(x) = 1 —— Qb(x) (2.13) 

UJ + 7T{X) 

is a solution of L c c(x) = 0. Similarly, if L c c(x) = 0, then 

b{x) = 1 —— Q ] c{x) (2.14) 

uj — 7r(x) 

solves L b b(x) = 0. 

Thus, in particular, given a pair {bi(x), b 2 (x)} of independent fundamental so- 
lutions of L b b(x) = 0, we can obtain from it a pair {ci(x), c 2 (x)} of independent 
fundamental solutions of L c c{x) = by using (|2.13|) , and vice versa. Therefore, with 
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no loss of generality, we henceforth assume, that the two pairs of independent funda- 
mental solutions {bi(x),b 2 (x)} and {c%(x), c 2 (x)}, are related by ( |2.13|) and (|2.14j) . 

The coefficients A^l(k), Af\(k) in ( |2.9[ ) are to be adjusted according to ( |2.13| ) and 
(|2.14|) , and this was the purpose of introducing them in the first place. 

Thus, with no loss of generality, we may make the standard choice 

A$> = = 1 (2-15) 

in ( |2.9D . The coefficients are then determined by ( |2.13j ): 

(1) _ a(-oo)-ik 

c — 7 \ 

7T{ — 00) + UJ 

4 2) = ff (°°\ + ik . (2.16) 

7r(oo) + UJ 

We note that these b(x) c(x) mappings can break only if 

Qb = or g t c = 0, (2.17) 



for b(x) or c(x) that solve ([2.12 ). Do such solutions exist? Let us assume, for example, 



that Qb = and that L^b = 0. From the first equation in (|2.12j ) (or in (|2.6|) ), we 
see that this is possible if and only if u ± ir(x) = 0, which clearly cannot hold if 
d x 7r(x) 7^ 0. A similar argument holds for c = 0. Thus, if d x 7t(x) ^ 0, the 
mappings ( |2.13j ) and (|2.14j ) are one-to-one. In particular, a bound state in Lb implies 
a bound state in L c (at the same energy) and vice-versa. 



An interesting related result concerns the wronskians and W c . From ( 2.11 



and from ( ^.13| ) and ( |2.14[ ) it follows immediately that for pairs of independent fun- 



damental solutions {bi(x) , b 2 (x)} and {ci(x), c 2 (x)} we have 

rrr c 2 d x ci - cid x c 2 b 2 d x b x - b x d x b 2 

W c = — = Ci& 2 - c 2 6i = ■ — p-r = W b . 2.18 

UJ — TX[X) UJ + TC(X) 

The wronskians of pairs of independent fundamental solutions of L b and L c , which 
are related via ( |2.13| ) and ( |2.14| ) are equal! 



To summarize, if d x 7r(x) ^ 0, L b and L c have the same set of energy eigenvalues 
and their eigenfunctions are in one-to-one correspondence. 
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If, however, ir =const., then we are back to the familiar "supersymmetric" factor- 
ization 

Q^Qb={u 2 -n 2 )b, QQ ] c= {uj 2 -7t 2 )c, (2.19) 

and mappings 

c(x) = — *— Q b(x) , b(x) = —^—Q ] c(x) . (2.20) 

UJ + TT OJ — It 

As is well known from the literature on supersymmetric quantum mechanics, the 
mappings (|2.20|) break down if either Qb = or Q^c = 0, in which case the two 
operators Q^Q and QQ^ are isospectral, but only up to a "zero-mode" (or rather, an 
oj 2 = tt 2 mode), which belongs to the spectrum of only one of the operators^. The 
case 7r(x) = brings us back to the GN model. Supersymmetric quantum mechanical 
considerations were quite useful in the study of fermion bags in [|Kj . 

The "Witten index" associated with the pair of isospectral operators Lf, and L c , is 
always null for backgrounds in which d x ir(x) ^ 0, since they are absolutely isospectral, 
and not only up to zero modes. There is no interesting topology associated with 
spectral mismatches of Lb and L c . This is not surprising at all, since the NJL model, 
with its continuous axial symmetry, does not support topological solitons. This is in 
contrast to the GN model, for which tt = 0, which contains topological kinks, whose 
topological charge is essentially the Witten index of the pair of operators ( |2.19| ). 

We note in passing that isospectrality of Lf, and L c which we have just proved, 
is consistent with the 75 symmetry of the system of equations in ( |2.4j ), which relates 
the resolvent of D with that of D — — 75.D75. Due to this symmetry, we can map the 
pair of equations L b b(x, y) = 5(x — y) and L c c(x, y) = —5(x — y) (Eqs. ([2.6|)) on each 
other by 

b(x,y) «-> —c(x,y) together with (cr, n) — > (—a, —n) . (2-21) 

(Note that under these reflections we also have a(x, y) «-> —d(x, y), as we can see from 
(|2.5|).) The reflection (cr, tt) — > (— cr, — n) just shifts both asymptotic chiral angles 9± 

3 This is true for short range decaying potentials on the whole real line. Strictly speaking, (to the 
best of our knowledge) only the case tt = appears in the literature on supersymmetric quantum 
mechanics. 
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by the same amount ir, and clearly does not change the physics. Since this reflection 
interchanges b(x, y) and c(x, y) without affecting the physics, these two objects must 
have the same singularities as functions of a>, consistent with isospectrality of L b and 



2.2 The Diagonal Resolvent 

Following [18|, |Tl| we define the diagonal resolvent (x \iD~ l \x) symmetrically as 

/ A{x) B(x) 

v C(x) D(x) 



(x\ — iD \x) 



(2.22) 



1 / a(x,y)+a(y,x) b(x,y) + b(y,x) 
= - lim 

V c(x, y) + c(y, x) d(x,y) + d(y,x) , q _ r , 

Here A(x) through D(x) stand for the entries of the diagonal resolvent, which follow- 
ing ( j2.5|) and (|2.8p have the compact representationF] 

bx(x)b 2 (x) 



B(x) 
C{x) 



ci(x)c 2 (x) 



D(x) 



A(x) 



1 % + 2a{x)} B (g) 

2 uj + it(x) 

1 [d x - 2a(x)} C (x) 

2 uj — ir(x) 



(2.23) 



We now use the generalized "supersymmetry" of the Dirac operator, which we 
discussed in the previous subsection, to deduce some important properties of the 
functions A(x) through D(x). 

From (|2.23| ) and from ( |2.3|) we we have 



A(x) = -- 



id x -2a(x){ c 1 c 2 \ i c 2 Q^ ci + CiQ 1 c 2 



2 lj-tt(x) V W c ) 2W C lj-tt(x) 



Using ( 2.14j) first, and then ( |2.13 ), we rewrite this expression as 



A(x) 



2W, 



(c 2 bx + C1&2) 



% b x Qb 2 + b 2 Qbi 
2W C uj + 7c(x) 



4 ^4, B, C and D are obviously functions of ui as well. For notational simplicity we suppress their 
explicit uj dependence. 
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Then, using the fact that W c = W b (Eq. (gTTg )) and (|2~23l ), we rewrite the last 
expression as 



A(x) 



i d x + 2a bib 2 \ i (d x + 2a)B{x) 



(2.24) 



2u + 7r(a;) \ W b J 2 u; + ir(x) 

Thus, finally, 

A{x) = D(x) . 

Super symmetry renders the diagonal elements A and D equal. 

Due to ( |2.23|) , A = D is also a first order differential equation relating B and C. 
We can also relate the off diagonal elements B and C to each other more directly. 
From (|2.23|) and from ( |2.13|) we find 

rW __^_ (Qh)(Qb 2 ) 

C[X} - W c ~ ( w + 7r) 2 W c - 1 } 

After some algebra, and using ( [2.18|) , we can rewrite this as 



-(« + nfC = a 2 B + aB' + |^ 

The combination b'^/Wb appears in B" = (b^/Wb)" '. After using Lbb\^ = to 
eliminate b'{ and b' 2 ' from S", we find 

(JIT 1 ^ 



6i«/ 2 1 „„ tt'B' 



-B" 



W b 2 2(w + tt) 
Thus, finally, we have 

- (u; + tt) 2 C = -B" , - B' 

1 ; 2 T \ 2{u + tt) J 

In a similar manner we can prove that 



2,2 / 2 , 
CT +7T — u — UO + 



Co> + 7T 



5 



^ _ a f _ ^ + 



an 



{uo - n) 2 B = — C" + a 



TT 



C" 



2 , / 2 
7T + CT — C<J 



LO + 7T 



(77T 



B. (2.26) 



C. (2.27) 



2 y 2{ui — 7r) y \ a; — 7r / 

We can simplify (|2.26|) and ( |2.27|) further. After some algebra, and using ( 2.23|) 
we arrive at 

C(x) = l - rT d x D(x) - U ~ B(x) 

B(x) = : 1 dx A(x)- UJ + n{ , X \ C(x) 



(2.28) 
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Supersymmetry, namely, isospectrality of L& and L c , enables us to relate the diagonal 
resolvents of these operators, B and C, to each other. 

Thus, we can use ( |2.23j ), ( |2.24| ) and ( |2.28| ) to eliminate three of the entries of the 



diagonal resolvent in ( |2.23|) , in terms of the fourth. 

Note that the two relations in (12.281) transform into each other under 



B «-> — C simultaneously with (er, it) 



-er, -7r 



(2.29) 



in consistency with ( 2.21|) . The relations in (|2.28 ) are linear and homogeneous, with 
coefficients that for d x n(x) ^ do not introduce additional singularities in the uj 
plane. Thus, we see, once more, that B and C have the same singularities in the uj 



plane. We refer the reader to Section 4 in for concrete examples of such resolvents. 

The case tt(x) = brings us back to the GN model. In the GN model, our B 
and C, coincide, respectively, with ujR- and —ujR + , defined in Eqs. (9) and (10) in 
10 1 . With these identifications, the relation A = D (Eq. ( j2.24|) ) coincides essentially 



with Eq. (18) of Jl0|. The relations (|2.26|) and ( [2.27|) were not discussed in [fL0~l , but 
one can verify them, for example, for the resolvents corresponding to the kink case 
cr(x) = mtanhmi (Eq. (29) in |l(|), for which 

2 



C 



UJ 



UJ 



2 ' 



B 



m sech mx 



- 1 



UJ 



C. 
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3 Bilinear Fermion Condensates and Vanishing of 
the Spatial Fermion Current 



Following basic principles of quantum field theory, we may write the most generic 
flavor- singlet bilinear fermion condensate in our static background as 



reg 



A , . dui 
N — tr 
2vr 



r(x| 



\x) 



(3.1) 



A(x) B{x) \ ( A B 

C(x) D(x) ) \CD 
\ i \ / / \ / VAC _ 

where we have used ( p. 22 ). Here a = 1, • • • , N is a flavor index, and the trace is taken 
over Dirac indices a, (3. As usual, we regularized this condensate by subtracting from 
it a short distance divergent piece embodied here by the diagonal resolvent 

( A B \ 1 / imcosO uj + msin# 

(x | — iD \x ) VAC - 



C D j 



2Vm 2 - u 2 



VAC 



\ —to + msin^ imcosO 

(3.2) 

of the Dirac operator in a vacuum configuration u VAC = mcos8 and vr VAC = msin#. 
In our convention for 7 matrices ( |2.1|) we have 

I A(x) B(x) \ ^ A(x)+D(x) 1 A(x)-D(x) l | . B{x) - C(x) | B(x) + C(x) 
K C{x) D{x) ) 2 21 2 2 

(3.3) 

An important condensate is the expectation value of the fermion current (j M (x)). 
In particular, consider its spatial component. In our static background (<j(x),ti(x)), 
it must, of course, vanish identically 



(j\x)) = 0. 



(3.4) 



Thus, substituting T = 7 1 in (|3.1| ) and using ( |3.3|) we find 



( 3 \x))=iN j^- [A{x)-D{x)\ 



(3.5) 
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But we have already proved that A(x) = D(x) in any static background (a(x),7r(x)) 
(Eq. (|2.24|) ). Thus, each frequency component of (j 1 ) vanishes separately, and (|3.4|) 
holds identically. It is remarkable that the generalized supersymmetry of the Dirac 
operator guarantees the consistency of any static (a(x),7i(x)) background. 

Expressions for other bilinear condensates may be derived in a similar manner 
(here we write the unsubtracted quantities). Thus, substituting T = 7 in (|3.1|) and 
using ( p.3|) , ( p.24|) and ( [2. 28j ), we find that the fermion density is 

(j°(x)) =%N — [B(x) - C(x)) = iN / — \^—A 1 J • (3-6) 

Similarly, the scalar and pseudoscalar condensates are 

$(x)i>(x)) =Nj^ [A(x) + D(x)} = 2Nj^- D{x) , (3.7) 

and 

{^{x)^{x)) = N / — [B(x) + C(x)} = N / 1 ; V ^ • ( 3 - 8 ) 
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